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We study the effects of moderately strong magnetic fields on the properties of color-fiavor locked 
color superconducting quark matter in the framework of the Nambu-Jona-Lasinio model. We find 
that the energy gaps, which describe the color superconducting pairing as well as the magnetiza- 
tion, are oscillating functions of the magnetic field. Also, we observe that the oscillations of the 
magnetization can be so strong that homogeneous quark matter becomes metastable for a range of 
parameters. We suggest that this points to the possibility of magnetic domains or other types of 
magnetic inhomogeneities in the quark cores of magnetars. 
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I. INTRODUCTION 

It has long been suggested that the superdense core 
of compact stars may be composed of deconfined quark 
matter However, calculations using effective mod- 
els for the stroiigly coupled limit of quantum chromody- 
namics (QCD) 0,3 j ^ '^sll as ab initio calculations per- 
formed in weak coupling d, |^ , predict that a phase tran- 
sition between normal and color superconducting (CSC) 
quark matter takes place at high baryon densities and 
sufficiently low temperatures (for early papers on color 
superconductivity see Ref . ^ , and for recent reviews see 
Ref. Q). Therefore, it is plausible that the inner region 
of compact stars may consist of color superconducting 
matter. 

Strong magnetic fields exist on the surface of compact 
stars, e.g., B < 10^^ G for ordinary neutron stars 
while for magnetars they can be as large as B ~ 10^^ G 
0. Also, recent studies have given support to the idea 
that old neutron stars could have previously been magne- 
tars with magnetic fields that decayed over many years 
[Tol [Tlj . This then suggests that magnetars could also 
have color superconducting cores. Therefore, a detailed 
study of the effects of strong magnetic fields on color su- 
perconductors may be very important to understand the 
physics of magnetars. 

According to the standard magnetar model by Dun- 
can and Thompson , the energy bursts of soft gamma 
ray repeaters (SGR's) are caused by solid fractures in- 
duced by strong magnetic fields in the crust of the star. 
Many properties of SGR's are well described using the 
hadronic magnetar model [13i | but there are still some 
remaining issues, such as the quasi-periodicity of active 
phases, which require alternative ways of understanding 
them |14J . 

The physical upper limit for the magnetic field in a 
gravitationally bound star is B ~ 10^^ G [1, [lB|, which 
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is obtained by comparing the energy of magnetic and 
gravitational fields. However, if quark stars are self- 
bound rather than gravitationally bound objects this up- 
per limit could be even higher. The typical energy scale 
defined by these strong fields is of the same order of mag- 
nitude as the color superconducting gap present in the 
quasiparticles energy spectrum. Thus, one expects that 
strong magnetic fields should affect the diquark pairing. 
The first step in this direction was presented in Ref. [l6| , 
where the effects of an extremely strong magnetic field 
on the pairing of a color-flavor locked (CFL) supercon- 
ductor ^17j were studied. The medium magnetic field in 
the CFL phase is a linear combination of a gluon and the 
usual vacuum magnetic fields This "rotated" mag- 
netic field is not subject to the Meissner effect and, thus, 
freely penetrates the system. 

In this paper we investigate the effects of a strong mag- 
netic field (namely, eB/fi^ < 1 where is the baryon 
chemical potential of quarks) on color-flavor locked su- 
perconducting quark matter. This is expected to model 
the conditions which exist at the inner core of magnetars, 
where the density of matter is up to ten times higher 
than the nuclear saturation density. For simplicity's 
sake, we assume that such matter is sufficiently dense so 
that effects of non- vanishing quark masses (including the 
strange quark mass) could be neglected. Also, using the 
property of enforced neutrality of the color-flavor locked 
phase [101 , we assume that no electrons are needed in the 
magnetized version of the phase. By applying the same 
arguments as in Ref. jl^ . we expect that this is justified 
as soon as all quark quasiparticles remain gapped. In our 
analysis, therefore, the chemical potentials of up, down 
and strange quarks are equal. Admittedly, we neglect 
possible non-zero "color" chemical potentials that might 
be needed to enforce the color Gauss law How- 
ever, such chemical potentials are expected to be small 
for the range of parameters studied in this paper, i.e., 
eB/fi^ < 1. 

In order to determine the thermodynamic properties 
of the magnetic CFL (mCFL) phase, we solve the cor- 
responding gap equations for a range of magnetic fields. 
As expected, the magnetization of quark matter displays 
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de Haas- van Alphen oscillations and it can be as large 
as the applied magnetic field for a wide range of param- 
eters. This is in contrast to the magnetization of dense 
hadronic matter, which is usually considered to be neg- 
ligible ^21j . The large magnetization in magnctar mod- 
els with (either normal or color superconducting) quark 
matter cores could be related to physical properties that 
are distinctively different from models with only hadronic 
matter. 

In the next section we discuss the most important fea- 
tures concerning mCFL color superconductors and also 
illustrate the derivation of the free energy and the cor- 
responding gap equations (details of our calculations are 
explicitly given in Appendix A and B). In the same sec- 
tion we discuss the regularization scheme employed in 
our calculations and, in Sec. IIIH we show our numer- 
ical results. In the last section we give our conclu- 
sions. Throughout the paper the units are h = c — 
/cb = 1 and 4- vectors are denoted by capital letters, 
e.g., — {kQ,k). Also, the standard Minkowski metric 
Vfj-f = diag(+, -, -, -) is used. 



II. THE FREE ENERGY OF THE MAGNETIC 
CFL PHASE 

In cold, dense quark matter, color superconductivity is 
expected to occur due to attractive interactions between 
quarks located on the Fermi surface. At very low temper- 
atures and extreme high densities one can safely neglect 
the effects of nonzero quark masses and, in this case, 
the CFL phase is expected to be the true ground state of 
quark matter p^j . This phase is characterized by the con- 
densation of quark Cooper pairs in the color-antitriplet, 
flavor-antitriplet representation, which breaks the initial 
symmetry group [SU{3)]c(d SU{S) l ® S't/(3)K(^ C/(l)i3 
U{1)a into the diagonal SU{3)ii+L+c subgroup (note 
that U{1)a is an approximate symmetry in dense mat- 
ter) . After taking into account that all the gluons become 
massive through the Higgs mechanism, one finds that ten 
Nambu-Goldstone bosons appear due to the breaking of 
global symmetries [12, 1131 ■ Following the same conven- 
tion as in the QCD vacuum, these bosons are the tt^, 
7r°, if^, if^, rj, and r]' mesons, and if, which is a 
massless field related to the baryon symmetry breaking. 

The CFL phase has some very unusual properties, for 
instance, it is a baryon superfluid but not an electro- 
magnetic superconductor [17| . In fact, when it comes 
to its electromagnetic properties a CFL superconduc- 
tor is more adequately described as an isotropic elec- 
tromagnetic insulator |19|. In conventional supercon- 
ductivity, the C/(l) gauge symmetry of electromagnetism 
is spontaneously broken by the condensation of electron 
Cooper pairs, which gives an effective mass to the pho- 
ton that screens external magnetic fields (electromag- 
netic Meissner effect). On the other hand, in color su- 
perconductors the initial [C/(l)]em symmetry, whose gen- 
erator coincides with one of the vector-like generators of 



SU{3)L<E)SU{3)ii, is not really broken but "rotated" 
The new group, denoted here as [J7(l)]em, corresponds to 
a massless linear combination of the vacuum photon and 
the eighth gluon. This is analogous to the electromag- 
netic [C/(l)]em symmetry that remains unbroken after the 
electroweak symmetry breaking. Using the same conven- 
tion for colors as in Ref . [l^, , it can be shown [13, HE] 
that the "rotated" electric Q-charges of the quarks are 
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in units of the Q-charge of the electron, which is e = 
ecos9 with 9 = aiccos{g / + as the mixing an- 

gle [S]. In color-flavor space, the associated charge op- 
erator Q is given by Q — (5/(g)lc— I/® Qc where 
Qc = -\s/V3 = diag(-l/3, -1/3, 2/3). 

The presence of an external magnetic field already re- 
duces the symmetry of the model because only s and 
d quarks have the same electric charges, i.e., SU{3)l ® 
SU{3)r -> SU{2)l ® SUi2)ii. Therefore, the pattern of 
symmetry breaking that defines the mCFL phase is given 
by [SU{3)]c ® SU{2)l SU{2)r ® U{1)b ® U{1)a ® 
t/(l)^ =^ 5'[/(2)fl+i+c, where C/(l)^ corresponds to 
an anomaly-free current formed by a linear combination 
of the s, d, and u axial currents [261]. Note that in this 
case there are only six Nambu-Goldstone bosons (i.e., 7r°, 
, K'^, r], r]' , and (p) due to symmetry breaking and all 
of them are neutral with respect to the Q-charge. As 
expected, the masses of the charged Nambu-Goldstone 
bosons of the CFL phase (i.e. tt^ and K^) have terms 
that are proportional to \/eB (see also Ref. [28^). 

It is convenient at this point to discuss the energy 
scales that appear in our study in more detail. The rel- 
evant mass scales in the problem are the quark chem- 
ical potential fi, the temperature T, and the magnetic 
length I — 1/^/eB. For very high, but still realis- 
tic densities, fi ~ 0.5 GeV. If the magnetar's interior 
fields are indeed as large as i? ^ 10^* G the corre- 
sponding magnetic energy scale approaches the QCD 
scale, VeB < 77 MeV. (Note the following relation: 
~ 7.69 X 10~^y/B/lG MeV.) The relevant dimen- 
sionless parameter in the problem is the ratio ei?//i^, 
which is smaller than one even for the largest possible 
fields. Thus, from a phenomenological point of view it 
would suffice to study only the regime of small eB / fi^ . 
However, our analysis extends also to larger fields that 
correspond to eB/^^ ~ 1. Note that at such extremely 
high fields, the magnetic length is comparable or even 
smaller than the average interquark distance. Then, 
all usual assumptions about the color superconductor's 
in-medium electromagnetic properties should be rean- 
alyzed. In fact, in this limit the most important ef- 
fect caused by the field might be the so-called magnetic 
catalysis of chiral symmetry breaking, as in the vacuum 
[H, US • The corresponding new ground state is char- 
acterized by different condensates. Its symmetry, though. 
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is similar to the symmetry in the mCFL phase, except 
for the baryon number symmetry that remains unbroken 

We use a massless three-flavor quark model with a local 
NJL-type interaction to describe the mCFL phase. The 
free energy density of this system composed of quarks 
in the presence of an external magnetic field H is a func- 
tional that depends on the gap functions, the chemical 
potential, the temperature, and the induced magnetic 
field B. These two fields have different physical meaning 
and need to be distinguished. For instance, in an ideal- 
ized model for the core of a color superconducting magne- 
tar, H corresponds to the field present in the outer layers 
of the star. Inside the core one has to take into account 
the magnetization of the medium and in this case B has 
to be used. We assume that both H and B are uniform 
fields that point in the f-direction. This approximation 
is valid as long as the fields do not change appreciably 
within the relevant length scales defined throughout the 
computation of the system's thermodynamic quantities 
such as the quark number density. This is indeed the 
case here because the characteristic length scale associ- 
ated with the problem is of the order of one fermi. 

The partition function of the system is given by 
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where Af is the normalization constant. The Lagrangian 
density reads 



£^^{i^ + eQ4 + /i7o)^ + ji-^Pvi^cKi^cP^i^), (2) 

where A^^ describes the gauge field of the in-medium elec- 
tromagnetism [U{l)]em, /? = ^/T is the inverse temper- 
ature, and V is the 3-volume. The quark spinor car- 
ries color a = {b,g,r) = (1,2,3) and flavor a = (s,(i, it) 
indices, and the charge-conjugate spinors are defined as 
= C^p'^ and ijjc = ip'^C, where C = 17^7° is the charge 
conjugation matrix. We consider only the pairing in the 
antisymmetric channels and, thus, {Pri)'^p — ^75 (-""^^(-apri 
(no sum over rf) and = 7o-P^7o. The index rj labels 
the pairing channels, i.e., 77 = 1, 2, and 3 correspond to 
ud, us, and sd pairings, respectively. 

For each channel we introduce a complex scalar field 
(pn, with expectation value A,,. The four-fermion inter- 
action is bosonized via a Hubbard-Stratonovich transfor- 
mation, which then gives a Yukawa-type interaction, 

^(^F,^c)(V5cP,V') y(V^cP,V') + ^(MV'c)-^. 

(3) 

In the following, we neglect diquark fluctuations and, 
thus, set (pn = A^. Moreover, using the standard Nambu- 
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and the gap matrix <f>+ ~ ^'i]=i ^v^v ^^'^ Lagrangian 
density can be rewritten as 

CiX) = - ^ ^ + h{X)S-\X)^{X) , (5) 



ri=l 



where 



(6) 



and the following shorthand notation is used: 
-1 = [i^ + eQ4{X) ± fijo] and $- = 7o($+)Wo. 



0(Q)' 

The mCFL pairing is characterized by the residual 
SU{2)ii+L+c global symmetry, which means that the 
corresponding gap matrix is invariant under simultane- 
ous flavor (1^2) and color (1 2) exchanges. This 
implies that A = Ai = A2 and = A3. The CFL gap 
structure is recovered when A — (f). 

Even though all diquark pairs are chargeless with re- 
spect to the rotated electromagnetism they can be com- 
posed of either quarks with opposite Q-charges or Q- 
neutral quarks (in the following, when discussing quark 
charges we always have in mind the Q-charges even if this 
is not explicitly emphasized). The gap function (f> has 
only contributions from pairs of neutral quarks, while 
A is formed by pairs of quarks with opposite charges 
and pairs of neutral quarks. In the presence of a strong 
magnetic field one would naively expect that pairs made 
of quarks with opposite charges have smaller coherence 
lengths in comparison to the pairs formed only by neu- 
tral quarks, i.e., A should be larger than (j). In the next 
section we show that this is the case in the limit of very 
large fields, but is not generally true. 

The Gibbs free energy density Q of the mCFL phase is 
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where 



2A2 

^=^ + ^-r(T,M,A,0,B) 



(7) 



(8) 



G G 

The last term is the one- loop contribution of quarks, i.e.. 



r{T,fi,A,(t),B) = ilndet5"\ 



(9) 



The free energy of the vacuum is given by J-Vac — 

-rvac = -r(o,o,o,o,o). 

In equilibrium Q is evaluated at its stationary point 
with respect to A, ^, and B and it describes all the sys- 
tem's thermodynamic properties. The CFL free energy. 
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^0, has to be recovered when H 0. The stationary 
point corresponds to the solutions of the equations 



G fdT\ 
B = H + inM, 



(10a) 

(10b) 
(10c) 



where M is the magnetization of the system, which is 
given by M = (9r/5P)| stationary At the stationary 
point, the gaps depend explicitly on the induced field 
B. This field incorporates the magnetic properties of the 
medium, described by the magnetization M. 

In this paper we consider the zero temperature case 
only, i.e., T — because typical temperatures of matter 
in stellar cores are much less than 1 MeV [ll| . According 
to Eq. (jA33[) in Appendix \X[ the one- loop quark contri- 
bution to the free energy reads 



r(0,Ai,A,<^,B) 



= 3P(0) + F(Ai) 
+ 4F(A), 



^(A2) 



(11) 



where A1/2 = ^(\/0^ -I- 8A^ ± 0). The neutral quark 
contributions are given by the three terms containing P- 
functions 



P(0) = trs^o K(0) + So-(0)] , 



where the trace is defined as tr^^ 



J i2¥Fi- ■ 



(12) 



and 



= \/(p T A*)^ + 0^ with p = \/p§ -I- p]_. On the 
other hand, the one-loop contributions that come from 
charged quarks are given in terms of the following func- 
tion: 



P(A) = trB[E+{A)+E^{A)] 



(13) 



The trace here is defined as a sum over the Landau lev- 
els with an integral over the longitudinal momentum, 
trB[. • ■] = 1^ Er=o I-oo dpal ■ ■] with a„ = 2 - 6no- 
The quasiparticle dispersion relation of charged quarks is 
E^[A] = T m)^ + A2, where es = y^pl + 2eBn. 

Using the well-known Euler-McLaurin summation for- 
mula 



oc 
n=0 



1 

dxfix)--f'{0) + --- , (14) 



we see that, formally, Huib^q F{x) — P{x). Therefore, 
after taking the limit A = (p and B — H — we recover 
the free energy of the CFL phase in the absence of an 
external field. 

The traces in the definition of the P and F functions in- 
volve integrations and sums over the whole phase space, 
which means that these functions diverge in the ultra- 
violet. In order to obtain the physically meaningful free 
energy density of the system, these functions must be reg- 
ularized. In effective quark models such as the Nambu- 
Jona-Lasinio model used here, it is common to restrict 



the phase space by introducing a finite sharp cutoff in 
momentum space. However, because of the special prop- 
erties of the system in a magnetic field such a prescription 
is not very useful. Utilizing a sharp cutoff when an en- 
ergy spectrum with discrete Landau levels is considered 
would introduce unphysical discontinuities in many ther- 
modynamical quantities. 

In this study, therefore, we regularize the traces in 
Eqs. (fT2l) and p3| by introducing a smooth cutoff func- 
tion h\ (where A is a constant with the dimension of 
energy). The cutoff function should approach 1 at small 
energies (i.e., £ ^ A) and at large energies (i.e., e ^ A). 
Providing that /ia falls off sufficiently fast in the ultravi- 
olet, the regularized functions, 



d?p 
(2^ 



and 



^a(A) 



eB ^ 



n=0 



hK[E+{^)+E^{4>)\. (15) 

O 

dp3/iA [E+{A)+Eg{A)\ . 

(16) 



are free from divergences. 

After using several different choices for /ia we decided 
on the Gaussian-like form, i.e., 



/lA = exp (-£^/A^) , 



(17) 



where e = p in Eq. (fT5|) and e = £s in Eq. (fT6| . Note 
that a sharp cutoff could be implemented with the func- 
tion /lA — (A — e) [where 9{x) is the step function]. 
While this would produce no apparent abnormalities in 
the behavior of Pa(0), it leads to unphysical discontinu- 
ities in Pa (A) and, consequently, in the free energy as a 
function of the magnetic field. 

In the numerical calculation in Sec. IIIH we use the 
cutoff function in Eq. ([17]) with A = 1 GeV. It should be 
mentioned here that we have checked the cutoff indepen- 
dence of our qualitative results by varying the value of 
A from 1 GcV to 2 GeV, and simultaneously readjusting 
the diquark coupling constant G so as to keep the CFL 
gap at i? = unchanged. To avoid a potential confusion, 
let us additionally mention that the dependence of the 
results on the cuttoff parameter A cannot be completely 
eliminated in the analysis. This is a general feature of all 
non-renormalizable effective models. 

When using the smooth cutoff function ([T7]) in our cal- 
culation, we get no divergences in the theory. Because no 
consistent continuum limit (A — > 00) exists in the model 
at hand, there is also no need for any additional (finite) 
renormalization and/or subtractions in our calculations. 
Strictly speaking, one should only remember that all 
physical parameters of the model (e.g., the charge and the 
magnetic field) are defined at the scale of the cutoff (i.e., 
1 GeV for our choice of cutoff A). Because of the slow 
logarithmic running in the complete gauge theory, the 
corresponding values of parameters differ by less then 5% 
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from the parameters defined at the conventional infrared 
scale set by the electron mass. Also, since most of the 
results of interest depend on the renormalization group 
invariant quantity eS, we choose to ignore these minor 
corrections. [Alternatively, one might consider an explicit 
finite multiplicative renormalization cr = e/vl + Ct? 
and Br = By/l + Ce^, with C cx \n{A'^ /ml) |31i]. Here 
the new quantities, cr and _Br, are defined at the scale 
set by the electron mass. However, we do not find this 
necessary for the purposes of the current study.] 

By making use of the regularized functions Pa{4') and 
Fa (A), the gap equations (jlOap and (jlObp become 



A 



G 

T 

G 



[i?iPi(Ai) + i?2Pi(A2) + 4AFi(A)] ,(18a) 



^ = -[30Pi(0) + C/iPi(Ai) + C/2Pi(A2)], (18b) 



where the functions Pi and Pi are defined as follows 

1 dP^{4>) 



Pi (A) 



^ 1 rfPA(A) 

~ A dA ■ 



19a) 



(19b) 



In the gap equations above we also introduced the nota- 
tion 



(20a) 
(20b) 



Pi/2 = 2A 1± 



Ul/2 



± 



8A' 



+ 8A2 



The gap equations derived here are slightly different than 
those obtained in Refs. [H, H^]. The reason is likely to 
be the difference in the structure of the interaction be- 
tween the quarks. However, as shown in the next section, 
at asymptotically high fields our numerical solutions for 
the gap functions are in qualitative agreement with the 
analytical formulas obtained in Refs. [l^, H^l • 



III. NUMERICAL RESULTS 

The free parameters of the model (G and A) were 
set to yield a CFL gap of either 0o — 10 MeV or 
(j}o = 25 MeV at M = 500 MeV when P = 0. (All nu- 
merical results below correspond to the following specific 
choice of the coupling constant and the cutoff parameter: 
G = 4.32 GeV"^ or G = 5.15 GeV"^ and A = 1 GeV.) 
In Fig. [1] we plotted the ratio between the mCFL gaps 
and the CFL gap as functions of eB/fi^. As was men- 
tioned in the last section, we have checked the robust- 
ness of our results by changing the parameters A and G. 
When eB/ fj? > 0.3 the gap A, which receives contribu- 
tions from pairs of chargeless quarks and also pairs with 
opposite charges, differs considerably from (/), where only 
chargeless quarks enter in the pairing. For smaller fields 



(eP//i2 < 0.1) the mCFL gaps are practically the same 
as 4>a. 

For ultrastrong fields {eB/jj? > 2) A is much larger 
than (f). This is consistent with the analytical solutions 
of the mCFL gap equations found in Refs. [iBjIl^l- There 
the authors considered fields so large that only the lowest 
Landau level contributes to the dynamics. As argued in 
the previous section, such fields should already probe the 
dynamics of the vacuum and the dynamics due to the 
magnetic catalysis should be taken into account [l^, [2^, 
Isoj . Moreover, the corresponding field strengths appear 
to be of the order of P ~ 8.5 x 10^^ G, assuming /i = 500 
MeV, which may be too large to be found in compact 
stars [32j . 

In Fig. [T] we see that the gaps display magnetic oscilla- 
tions with respect to eB/fj,^. These oscillations share the 
same physical origin as the de Haas- van Alphen oscilla- 
tions observed in metals. They appear as a consequence 
of the oscillatory structure in the density of states, which 
is imposed by the quantization of the energy levels asso- 
ciated the orbital motion of charged particles [s^ . 

At zero temperature most electronic properties of met- 
als depend on the density of states on the Fermi surface. 
Therefore an oscillatory behaviour as a function of B 
should appear in any quantity that depends on the den- 
sity of states on the Fermi energy. In the case of the 
mCFL color superconductor, it is the quark density of 
states on the Fermi surface that is relevant. Since the 
gaps in the excitation spectrum depend on the quark 
density of states and every physical quantity that we are 
interested in depends on the gaps, some type of mag- 
netic oscillations should appear. The presence of nonzero 
gaps, however, smears out the Fermi surface so that the 
oscillatory structure is considerably reduced. This then 
explains the smoothness of the oscillations and their de- 
pendence on the magnitude of the gaps in Fig. [TJ 
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FIG. 1: (Color online) Ratios A/4>o and 0/0o versus eB / 
for two sets of parameters that yield 4>o = 10 MeV and (po = 
25 MeV. 

In the presence of strong magnetic fields, the magneti- 
zation M would affect the properties of the mCFL phase. 
As mentioned in the introduction, this might be relevant 
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to explain the magnetic properties of magnetars with 
color superconducting quark cores. A phenomenon that 
is quite often observed in magnetic systems is the for- 
mation of magnetic domains ^33]. As in metals, the large 
magnitude de Haas- van Alphen oscillations in the magne- 
tization can create regions where (dH/dB)^ < 0. These 
correspond to unstable or metastable states. The condi- 
tion for thermodynamic stability {dH/dB)^^ > implies 
that 47r |x(At, B)\ < 1, where x(m, B) = {dM/dB)^ is the 
differential susceptibility. When the differential suscep- 
tibility exceeds l/47r, which depends on the geometry of 
the system, a transition into a magnetic domain configu- 
ration may occur. The presence of magnetic domains in 
the crust of a neutron star was discussed by Blandford 
and Hernquist in Ref. [s^ . 

In Ref. [2l| it was shown that the magnetization of 
hadronic matter is negligible even for magnetar condi- 
tions, i.e., A-kM/B < 1 for B < lO^^G. In the case of 
color superconducting quark matter, however, the situa- 
tion is very different. In Fig.[2]we plotted this ratio versus 
ei3//i^ for a mCFL superconductor. The magnetization 
is significantly larger in this case and it displays de Haas- 
van Alphen oscillations that have a very large magnitude. 
Note that nonzero energy gaps make the Fermi surface 
look fuzzy and the magnetization's oscillations in mCFL 
quark matter are much smoother than those shown by 
the magnetization of cold, unpaired quark matter (35| . 
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FIG. 2: Ratio A-k M/B versus eB/^^ for two sets of parame- 
ters that yield 4>o = 10 MeV and (t)o = 25 MeV. 

Our results for the H-B curve are shown in Fig. [31 Sev- 
eral regions of thermodynamic instability are obtained for 
eB/iJ? < I. The metastable regions in Fig. [31 which cor- 
respond to (dH/dB)^ < 0, can be filtered out by either 
using a Maxwell construction or a mixed phase where mi- 
croscopic domains with nonequal magnetizations coexist. 
This could lead to several physical possibilities. As the 
field H increases the system could undergo successive 
phase transitions with discontinuous changes of the in- 
duced magnetic field B [s^l . When the mixed phases are 
formed, the relative size of domains with different mag- 
netizations would change with H so as to keep the aver- 
age induced magnetic field B continuous. In either case. 



since the magnitude of the fields involved is enormous, 
the system could potentially release an immense amount 
of energy. Further studies are needed to see whether this 
finding could have any important implications for mag- 
netars. 




0.1 0.150.2 0.3 0.5 0.7 1 1.5 2 



FIG. 3: External field eH/fj,^ versus eB/^^ for two sets of 
parameters that yield (j)o = 10 MeV and (po = 25 MeV. 



IV. CONCLUSIONS 

In this paper we studied the effects of a moderately 
strong magnetic field on the Cooper pairing dynamics in 
cold and dense three-fiavor quark matter. We solved the 
corresponding gap equations and calculated the magneti- 
zation for a wide range of magnetic fields, eJB//Lt^ < 1. We 
found that, as the magnetic field increases, the system 
undergoes a continuous crossover from the usual CFL 
phase to the mCFL phase. Notably, for eB/fi^ < 0.1, 
which corresponds to B < 4.2 x 10^^ G (provided that 
/i ~ 500 MeV), there is no large difference between the 
mCFL and CFL gaps. 

In this study we also showed that the gaps display 
magnetic oscillations, which is a direct consequence of 
the Landau quantization of the energy levels in a mag- 
netic field. Similar magnetic oscillations were predicted 
and later observed in type-II electronic superconductors 
[3^ . The effects these oscillations have on the transport 
properties of mCFL superconductors still remain to be 
understood. 

We showed that the magnetization of mCFL quark 
matter displays de Haas-van Alphen oscillations whose 
amplitude can be as large as the magnetic field for a 
wide range of magnetic fields. Since the magnetization 
of dense hadronic matter is negligible, the nonzero mag- 
netization of color superconducting quark matter may 
provide new observable effects that can help to distin- 
guish purely hadronic magnetars from color supercon- 
ducting ones. Our results for the oscillations of the mag- 
netization suggest that homogeneous quark matter may 
become metastable for the range of parameters that are 
phenomenologically relevant to magnetars. Because of 
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this there is a possibihty that magnetic domains or other 
magnetic inhomogeneities can be found in the quark cores 
of magnetars. The successive phase transitions com- 
ing from discontinuous changes of the induced magnetic 
field B during the stellar evolution can release a vast 
amount of energy that would heat up the star, which 
would then cool down by for example the emission of 
neutrinos. Therefore, bursts of neutrinos coming from 
magnetars with color superconducting cores even after 
the deleptonization period could be expected. 

It has been recently pointed out in Ref. [s^l that 
the chromomagnetic instability [ssj . which is present in 
neutral two-flavor color superconductivity at moderate 
densities [s^, can be removed by the formation of an 
inhomogeneous condensate of charged gluons and the 
corresponding induction of a magnetic field. The field 
strengths used in their approach are comparable to those 
that cause the large magnitude de Haas-van Alphen os- 
cillations. Both effects and/or their competition may be 
relevant to the explanation of strong magnetic fields and 
other unusual observed properties of magnetars. 

There are several questions regarding color supercon- 
ducting matter in strong fields and its consequent use in 
the study of superconducting magnetars that remain to 
be answered. One could, for example, study the cooling 
of their ultramagnetized cores. Also, the effects of strong 
magnetic fields on the global structure of the star, such 
as its mass and radius, are also very important and will 
be discussed elsewhere 0|. The consideration of the ef- 
fects of a nonzero strange quark mass on the low-energy 
description of mCFL color superconductors should also 
be investigated. 

The low-energy effective theory for the Nambu- 
Goldstone bosons in the mCFL phase was recently de- 
rived by Ferrer and Incera in Ref. 27[ , although the pion 
decay constant and the meson maximum velocities still 
have to be determined within the microscopic theory. Us- 
ing the zero field values of these parameters calculated in 
Ref. [1^ they obtained that the charged mesons decou- 
ple from the low-energy theory only when eB ^ 12 0q. 
Incidentally, this corresponds to the region where the os- 
cillations of the magnetization and the gaps become al- 
ready noticeable. This brings up an interesting question 
regarding the possible effects of the oscillations on the 
low-energy dynamics in the mCFL phase. This problem 
is left for future study. 

Note added. While finishing our paper we learned 
that a partially overlapping study was done by Kenji 
Fukushima and Harmen J. Warringa [41]. 
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APPENDIX A: EVALUATION OF THE 
DETERMINANT 

In this appendix we present the details of the computa- 
tion of r(T, /z, A, (j>, B). The usual way to compute these 
determinants is to use the identity 



det I ^ ^ j = det {AD - ACA'^B) (Al) 



for the determinant of a block matrix. Using this identity 
we can rewrite Eq. ^ as 



r = i In det ( [6+- J 

^ 0(Q)J ^ 0(Q)J J ' ^ ^ 



0(-Q)J 



where [G 



± 1-1 



0(Q) 
ll CO 

cates the calculations 



eQj^{X) ± /i7o]. However, the 



nontrivial color- flavor structure of \G^, ~J ^ compli- 

^ 0(±Q)J ^ 



1. Introducing the charge projectors 

The easiest way to calculate Eq. (|A2|) is to go back to 
the Lagrangian density and introduce the charge projec- 
tors in color-flavor space 



ri(o) — 

^(+) = 



diag(l,l, 0,1,1,0,0,0,1), (A3a) 
diag(0, 0,0, 0,0, 0,1,1,0), (A3b) 
diag(0,0,l,0,0,l,0,0,0), (A3c) 



which satisfy 



f^(a)f^(6) = <5abO(t), a, 5 =0,-1-,-, ^ r2(a) = 1. 

a=0,± 

(A4) 

In terms of these projectors the charge operator in color- 
flavor space la Q — — However, it is convenient 
to define the following charge operator in Nambu-Gorkov 
space Qng = diag^g(Q, —Q), with eigenstates 



where 




(A5) 



= , , V^(a) = ^(a)^- (A6) 

Wc{-a) J 

The charge operators can also be generalized to include 
the Nambu-Gorkov structure present throughout this 
derivation. In fact, one has 



NG 
(a) 



V 



(A7) 
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which then satisfies 

a=0,± 

(A8) 

where 'i'(a) = ^faf ^- Using these projectors Eq. ^ can 
be expressed as 



^~ ^ G ^ , 2 



i)=l a,b=0,± 

However, it can be shown that 



(A9) 



^l^S-^^^^ = 5.,17f,f 5-if^f„f , (AlO) 

where we used that 51a$+ri_a is the only combina- 
tion that is not identically zero. We define S'Q^ = 

VL^9 S^^{X\ which then implies that 



(All) 



a=0,± 



It is possible to show that ,S^^^\ = 0. 

The term involving the full inverse propagator in the 
Lagrangian density can be written as 

a=0,± 

= I E ^ia)S{a]^ia)- (A12) 



a=0,± 



Therefore, the Lagrangian density in Eq. ([5]) can be 
rewritten as 



where 



C / ^ 2 

Tj=l a=0,± 



^-1 ^ (i'^tia)] ^ *(a) 



(A14) 



[^o(a)] ^ = ± Ai7o, and = + ae^. The new 
gap matrices in the equation above are given by ) = 
$+rj(a) whereas = 7o(*fa))^7o- 



2. Lagrangian density in momentum space 

It is very convenient to express the Lagrangian den- 
sity in Eq. (jA9p in momentum space because it simpli- 



fies the c omp utation of the determinants. Here we fol- 
low Ref. \24 and use the method originally developed 
for charged fermions by Ritus in Ref. [431. -R-^f- Hli 
the method was also extended to include charged vector 
fields. In this approach the diagonalization of the Green's 
functions of charged fermions in a uniform magnetic field 
in momentum space is obtained using the eigenfunction 
matrices Ep{X). These functions are the wave functions 
corresponding to asymptotic states of charged fermions 
in a uniform background magnetic field. 

The EJX) functions were described in detail in 
Refs. [23.143I| and some of their main properties are pre- 
sented in Appendix B. Using these eigenfunctions it is 
possible to express the charged fields V'(±) as 

^(±){X) = E^^H^Mi)!^)), (A15a) 
i'(±){X) = E'^(±)(%))^fH^), (A15b) 

where, by definition, e'^^\x) = 70 (i^jl^^^ 70, 
P{±) = (Po, 0, ±V2eBn,p3), and n = 0, 1, 2, . . ., denotes 
the Landau levels. The 4-vector potential is in the Lan- 
dau gauge, i.e., A^^ = (0,0, Bx,0). Moreover, one can 
show that 

[G±^)]-i(X)4±)(X) = (±) ±^7o]. (A16) 

This implies that 



0(a) 



(a) 



\G, 



0(a) 



(a) _ 



*(a)(Aa) 



(A17) 



for a = 0,±. Also, we have defined P(q) = {pq,p) and 
[^o(o)]"^(Ao)) [^(0) ± f^lo]- In the following section 
we will compute the determinants using the momentum 
representation of the inverse propagators. 



3. Calculating the determinants 

The operators S^^ and S^^J^ are defined in a 72- 
dimensional vector space E = Ccj^'Djjirac^J^NG- How- 
ever, the charge projectors express E as the direct sum of 
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three different spaces, i.e., S = S(o) ® S(+) ® with 
vector bases given by 



-(0) 



{|si,±), |S2,±), |di,±), \d2,±), \u3,±)}, (A18) 



This matrix can be easily diagonahzed and its 
eigenvalues are Af , Aj, where Af 2 = 

i(yi^pT8|A|2 ± |(/,|)2. Moreover, the other inverse 
propagators are given by 



S+ => {!«!,+), |U2, +), |S3, -), Ms, -)}, (A19) 



S- ^ {|ui,-),|u2,-),|s3,+),|d3,+)}, 



(A20) 



where {|+),|— )} is the basis of Nambu-Gorkov space. 
The color, flavor, and Dirac structures are automatically 
taken into account by describing a quark spinor -0 as 



'^(±) 



where 



2x2 



(±) 



(±) 



o(±)J 



, (A25) 



1-2x2 , 



4' = {si,S2,S3,di,d2,d3,Ui,U2,U3}, 



(A21) 



where (1,2,3) = {b,g,r) denotes the color indices. Thus, 
S(o), and are vector spaces with dimensions 

40, 16, and 16, respectively. 

We can now compute the determinant of the inverse 
propagator in Eq. ^ in terms of its corresponding charge 
projections. In fact, one sees that the determinant splits 
into three separate pieces 

r(T,/z,A,0,i?)= r(a)= E ^Indet^^-^. 

0=0, ± a=0,± 

(A22) 

Note that in the evaluation of the determinants only the 
projection of Sf^J^ on the corresponding lower dimensional 

subspaces ^{a), which we call 5^"), is relevant. The in- 
verse propagator that appears in the evaluation of 
the determinant is defined as 



whereas 



1^0(0)] 



1-5x5 



(0) 



[G, 



0(0)J 



(0) 

1(8)1 



(A23) 



5x5 , 



*(0) = il5 



/O (/) A\ 



0-0 

A 



A 0/ 



*(±) = -«75Al2x2. (A26) 
Now, each term in Eq. (IA22P can be simplified using 
the identity in Eq. (|A1|) . which then gives 



(a) 



■ In det 



[GoV)]"MGo(,)]-^+$(a)$(.) , (A27) 



where we used that [G^^)]-^*^^) = -^a^G^ia)]'' ■ 
Since only the absolute square of the gaps appears in 
the expressions, from now on we take A, > 0. 

In order to compute the determinants it is convenient 
to use the chiral and energy projectors defined in Ap- 
pendix B. First, it can be shown that 



K^r'^G^iar' = 5:[pg-(c|p(„)|±/.)^]A^,), (A28) 

c=± 



(A24) which then gives 



1 , f pI - i\P{a}\ - f^? - >^fa) 

(a) = ^^T^cja 2^ In I 

PO,P(a) V 



rp2 



PO.Pa 



'Po-(|P(a)l+Ai)'-A2 



(a) 



^2 



(A29) 



where —^fa) eigenvalues of the matrices ^^a)^*a)- 

Also, tr^j^^ is the remaining trace over the color, flavor, 
and chiral indices. Moreover, |p(o)| = \/pl +p1 and 

|p(-l-)| — \/p1 + 2eBn, where n > is an integer that la- 
bels the Landau levels, andpj^ is the momentum perpen- 
dicular to the field. Also, it is clear now that r(+) — r(_) . 



The sums in the equation above are defined as follows 



E /bo,P(o))=TE / ^fiiu^k.P). (A30) 



PO,P(0) 
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and 



with eigenfunctions 



P0,P(±) k n=0 

(A31) 

where / is an arbitrary function, and a — 2 — dno stands 
for the fact that Landau levels with n > are doubly 
degenerate. Since these sums do not converge some sort 
of cutoff procedure has to be used. The Matsubara sum 
can be evaluated using the identity 



Dii±)[Pi±)], (B3) 



T^ln 



rp2 



= |z| + 2Tln(l + e-l"l/^). (A32) 



Therefore, we can now write the one-loop quark con- 
tribution to the pressure of the mCFL phase as 



r = 3Pt(0) + ^^(Ai) + Pt(A2) + AFt{A). 
The function Pt{4>) is defined as follows 

Pt{c^) = trB=o[E+i(l>)+Eo {((>)] 



(A33) 



+2TtiB= 



ln(l + e-^'o/T) + + g-^'o/T) , (A34) 



where trB=o [. .] ^ / ■ E^ld^] = ^(^^TITT^^, 

with p = \/ pii + pj_ = \p\- The zero temperature terms 
are divergent and they are regularized in Sec. [ill The 
other function in Eq. (|A33p reads 

Ft (A) - trs [E+ (A) + Eg {A)] 



+2TtiB 



ln(l 



-E+/T 



)+ln(l + e-^«/^) (A35) 



where trs [. . .] = Y.n=o iZo ^P: 



], E^[A] 



\/{£b Tm)^ + with SB = \Jp\ + 2eBn. Note that 
limB_>o ^T(a^) — Pt{x) and, once the limit A — (\) and 
S = iJ = is taken, we recover the free energy of the 
CFL phase in the absence of an external field. 



APPENDIX B: PROPERTIES OF RITUS' 
EIGENFUNCTIONS 

The transformation functions i?p for positively (+) 
and negatively (-) charged fermions are obtained as so- 
lutions of the field dependent eigenvalue equation 

(n(±) • ^)E%{X) - i?|(X)(P(±) • 7), (Bl) 
where P{±) — {po, 0, ±^/2eBn,p3) and 

E$iX) = J2E$^iX)A[a], (B2) 



where Di{±)[p(±)] are parabolic cylinder functions with 
argument p(^ j.-^ defined as 



= V2eB{xi +p2/eB), 
and index ^ (±) given by 



l(±) 



n ± — — - , 
2 2 



;(±) = o,i,2. 



(B4) 



(B5) 



whereas n = 0,l,2,..., denotes the Landau levels and a 
is the spin projection that can take values ±1. Moreover, 
the normalization constant is 



Cz(±) = (4^ei3)i/Vv/U±y!- 



(B6) 



The spin matrices A[a] in Eq. (|B2p are spin projectors. 
They are defined as 



A[a] = diag((5^i,(5^_ 1,(5^1, 
and satisfy the following relations 



±1, 



(B7) 



A[±]t = A[±], 
A[±]A[±] = A[±], 
7||A[±] = A[±]7||, 



Af 



= 1, (B8a) 
A[=f]A[±]^0, (B8b) 
7iA[±] = A[t]7±-(B8c) 



In the equations above we used that 7|| = (7", 7^) and 

Massless quarks can be fully characterized by the chiral 
and energy projection operators. 



Pbx = 



A 



1±75 



1 ± 70 7 • P(a) 



(a) 



(B9) 

with a = 0,±. (BIO) 



respectively. Similarly to the free case, in the presence of 
a uniform magnetic field these two operators commute. 

In Ref. 4JI it was shown that only massless quarks with 
the same chirality pair in the spin zero channel. One can 
show that there are four different gap functions that de- 
scribe the possible pairings. However, in an NJL theory 
with the gap functions independent of the 3- momentum, 
the total number of independent gaps is reduced to one. 
This means that the gap for quasi-particles and quasi- 
antiparticles are the same, in spite of the chirality of the 
particles. Moreover, the Dirac structure of the gap ma- 
trix is simply given by C75, where C = 17^7'^ is the 
charge conjugation matrix. 
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